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Abstract. Consider an n x n Hermitean matrix valued stochastic process 
{Ht}t>o where the elements evolve according to Ornstein-Uhlenbeck processes. 
It is well known that the eigenvalues perform a so called Dyson Brownian mo- 
tion, that is they behave as Ornstein-Uhlenbeck processes conditioned never to 
intersect. 

In this paper we study not only the eigenvalues of the full matrix, but also 
the eigenvalues of all the principal minors. That is, the eigenvalues of the k x k 
■ in the upper left corner of Ht- If you project this process to a space-like path it 

is a determinantal process and we compute the kernel. This kernel contains the 
well known CUE minor kernel, |JN061 IOR06| and the Dyson Brownian motion 
kernel |NF98| as special cases. 

In the bulk scaling limit of this kernel it is possible to recover a time-dependent 
generalisation of Boutillier's bead kernel jBou09| . 

We also compute the kernel for a process of intertwined Brownian motions 
introduced by Warren in |War07] . That too is a determinantal process along 
spacelike paths. 



1. Introduction 

In a classic paper Dyson |Dys62| introduced a dynamics on random Hermitean 
matrices where each free matrix element evolves independently of all others. They 
each form an Ornstein-Uhlenbeck process, that is a Brownian motion with a drift 
toward zero. He successfully analysed the associated dynamics of the eigenvalues. 
Even the fact that the eigenvalues form a Markov process is highly non-trivial. 

This has become one of the most well studied models of random matrix the- 
ory. It is beyond the scope of this paper to completely survey the literature but 
some notable results are these. It has been found to be a limit of discrete random 
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walks conditioned never to intersect, for which correlation kernels have been found; 
see [NF981 [Joh02l IJoh03[ iJohOSbl IJohOSal lEKOSl IKT02[ IKT03j . Other work in- 
cludes |Ben08l [SpoSTj ITW04[ R^^MOS] . In particular, in |TW04[ partial 



differential equations were derived for the Dyson process and related processes. 

Actually there are two processes that have been called Dyson Brownian motion. 
Let a Hermitean matrix Bt evolve according to the transition density, for s < t, 

(1) PT{Bt e dB\Bs = B) = Ce-^^^^-''*-^)'/^!-''?-) 

where qts = e~^^~^^ and C is the normalisation constant that makes this a proba- 
bility density. Then this is a stationary process and its stationary measure is called 
the Gaussian Unitary Ensemble (GUE), see |Meh04| chapter 9]. For annxn GUE 
matrix B, 

(2) Pr(5 G dB) = Ce-^'^' dB. 

What this boils down to is that the elements on the diagonal are independent 
Gaussians with mean and variance |. For the off-diagonal elements the real 
and imaginary parts are independent Gaussians with mean and variance \. For 
an n X n Hermitean matrix B let eig-B = (Ai < • ■ • < A„) denote the vector of 
eigenvalues of B. The distribution of the eigenvalues is exactly 

n 

(3) Pr(eig B edX) = C JJ (A^ - Xjf JJ e"^? dK 

l<i<j<n i=l 

where C is a normalisation constant, see |Meh04t chapter 3]. The transition density 

for the eigenvalues is then, by the Harish-Chandra formula, 

(4) 

Pr(eigi?, G rfA| eigi?. = A) = det[e-(^--''-^'^^)^/(^— H dx,. 

^ ' i=l 

Here, A denotes the Vandermonde determinant. See |Joh05b] for a readable 
overview. This expression is a Doob /i-transform of a Karlin-McGregor determi- 
nant, see |Doo01trKM59j . That means that this can be interpreted probabilistically 
as n Ornstein-Uhlenbeck processes evolving in time conditioned never to intersect. 

Of course it is very natural to consider n pure Brownian motions (xi(t), . . . , 
x„(t))teR conditioned never to intersect. The transition density for that process 
would be 

(5) Pr(x(t) e dx I x{s) = x) = C^py dei[e-^''^-^^^^ JJ dxi. 

^''^^ i=i 

This process has also been called Dyson's Brownian motion, but is not stationary. 
It can be realised by a Hermitean matrix where the elements evolve as Brownian 
motions. 
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For the purpose of this article consider the following model which we shall call 
the Dyson Brownian minor process or (DBM process). Let (-Bt)feR+ be an x 
Hermitean matrix-valued stochastic process started at t = with Bq given by the 
GUE distribution. Let the process evolve with transition density given by ([T]). For 
n = 1, . . . , let -Bj^"^ be the n x n sub matrix in the upper left corner (principal 
minor) of Bt- We are interested in all the (^^^) eigenvalues of -Sf"^ for n = 1, . . . , 
A^. 

If A is an eigenvalue of sf"^ then we shall say that there is a particle at (n, t, A). 
In this way we can think of the DBM process as a point process, that is a measure 
on configurations of points or particles on the space {1, . . . , A^} x M x M. It turns 
out that this is, in the terminology of |BF08j , a determinantal point process along 
space-like paths. More precisely, that means the following. 

For notation we shall write that 

{true if n > n' 
true if n = n' and t < t' 
false otherwise. 

and 

in,t)>{n',t') := -n{{n,t) < {n' ,t')). 

Theorem 1.1. Take a sequence {(nj, Xj, tj)}jL]^ of levels, positions and times. Let 
them follow a space like path, which means that 

(6) < ti < t2 < ■ ■ ■ < 4, 

(7) rii > n2 > ■ ■ ■ > Uk. 

Then the density of the event that there is a particle at time ti on level ui at 
position Xi in the Dyson Brownian minor process is 

(8) p((?2i,xi,ti), . . . , {nk,Xk,tk)) = det[K^'^^{{ni,Xi,ti), inj,Xj,tj))]ij^-^ 
where 

(9) ir^«^((n,x,t),(n',x',t')) = 
' E ^^^^e''^*'"*^^:+K^)^:'-i-^(^')e-(^'^ for {n,t) > {n',t'), 

- E ^^^^e"'^*'-*^^:+K^)^:'+K^')e-(^'^ for {n,t) < {n',t'). 
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Furthermore 

(10) ir^^^((n,x,t),(n',x',0) = -(l)^'^^{{n,x,t),{n\x\t')) 



Off r,,n' „—u'^+2ux+v^—2vx' 

+ 2("'-«)/2^A^ du dv- ^ 



{2TTiy Jy e~(*'~*)t; — u 

where 

(11) 0^^^((n,x,t),(n',x',t')) = 

2(n-n')/2gn'(t'-t) /" ^n-n'^^ _ y)p,*,_,(y, x') rfy, z/ (n, t) < {u' , t'), 

Jr 

0, otherwise. 

The contours of integration are such that 7 encloses the pole at the origin and T 
goes from —ioo to ioo in such a way that \u\ < \v\ always, see Figured 

Here /i* , for n = 0, 1, . . . , is the normalised Hermite polynomial of order n, see 
Section Em /i* for n = —1, —2, . . . , are defined to be zero, p* is the transition 
density of an Ornstein-Uhlenbeck process, see (15U]1 . if" is the nth anti-derivative 
of the Dirac delta function, see (1471) . This Theorem will be proved in Section [7. 11[ 

The term space-like path for a path in space-time satisfying ([6]) and ([7j) was 
coined in |BF08j and, while the reason for using that name is not made clear, the 
terminology has become standard. 

It is quite clear that this process, along a space-like path, is a Markov process. 
It is known from Dyson |Dys62| that this process is Markovian on a fix level, 
that is for constant n. Baryshnikov |Bar01j observed that it is a Markov process 



for fixed time going down one level (say from n to n — 1). A space-like path is a 
combination of steps in time and steps going down one level, i.e. it is a combination 
of consecutive Markov steps. 

One beautiful construction which is due to Warren | War 07] is the following. 
Start a 1-dimensional Brownian motion, say {x{t))t>o at the origin at time t = 0. 
Then start two new processes, say (|/i(t))t>o and {y2(t))t>o respectively, one above 
and one below, respectively, the old one. They evolve as Brownian motions except 
that they are pushed up and down, respectively, by x. For details see |War07j . 
It is then a Theorem that yi and 1/2 together form a Dyson Brownian motion in 
the sense that their transition density is of the same form as The process 
can be continued: one can start three processes above, between and below yi and 
?/2- These three will then be a Dyson Brownian motion of three particles and so 
on. This process occurs as a scaling limit in the study of a certain random tiling 
model [Nor 10] . To reduce the amount of numerical factors floating around we shall 
in this paper consider the Warren process to be driven by Brownian motions with 
variance t/2 instead of standard Brownian motion. This is just a rescaling that is 
not important. 

For this model we can show a result analogous to that of Theorem 11.11 
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Theorem 1.2. Take a sequence {{ni,ti,Xi)}f^^ of times as in the previous Propo- 
sition. Then the density of the event that there is a particle at time ti on level Ui 
at position Xi in Warren's process is 

(12) p((^i, ti, xi), . . . , (nfc, tfc, Xk)) = det[K^{{ni, ti, Xi), {uj, tj, Xj))]ij^^ 
where 

(13) K^{{n,x,t),{n\x\t')) = 

71 E \/^^T7yr (0'^"+'(^/v^)^n'+/(^Vv^)e"(^')'/*', for in,t) > {n',t'), 

y "Tf E \/ ^^^"(^^:+/(^/v^)^:'+K^Vv^)e-(^')'/^ forM < {n',t'). 
Furthermore 

(14) 2^("-"')ir^((n, X, t), (n', x\ t')) = -0^((n, x, t), (n', x', t')) 

2 ^n/2ur\n' /2 r /• ^n' ^-u'^+2ux/ y^+v^ -2vx' /VF 

+ -1= du dv- 



(27ri)2 y/t 7p m"- tv — u 

where 
(15) 

lH"-\x-y)p,^,{y,x')dy, ^f {n,t) < {n',t'), 

0, otherwise. 

The contours of integration are such that 7 encloses the pole at the origin and T 
goes from —zoo to ioo in such a way that \u\ < \v\ always, see FigureUl 



Again /i* is the normahsed Hermite polynomial of order n, see Section 16.21 and 
if" is the nth anti-derivative of the Dirac delta function, see (1471) . p is the tran- 
sition density of a Brownian motion, see (l53l) . This Theorem will be proved in 
Section I7.11[ 

Note that although the kernels and K^^-^ are just a change of variables 
from each other, the underlying processes are different in an essential way. As 
noted in |ANvM10] , the difference of the eigenvalues of successive levels are pushed 
apart by a constant drift when they are close. In the construction due to Warren 
in |War07] the difference of the particles on successive levels behave like the abso- 
lute value of a Brownian motion when they are close. It is to be remarked that if 
the minor process was constructed with Brownian motions replacing the Ornstein- 
Uhlenbeck prcesses, the marginals along space-like paths would agree with the 
Warren process. 

No article on random matrices is complete without a scaling limit, so let us do 
one of those. In |Bou09j . Boutillier introduced a one-parameter family of models 
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which are point processes on N x M. On each individual copy of M, it speciahses to 
a determinantal process with the sine kernel which is so prevalent in all branches of 
random matrix theory, see |Meh04] . Furthermore, on successive lines the particles 
interlace. By that I mean that if there are particles at {n,xi) and (n, X2), then 
there exists a particle (n + 1, such that xi < y < X2 almost surely. As a scaling 
limit of K^^^ above we recover a kernel which specialises to the Boutillier Bead 
kernel at a fixed time. One way to interpret this is to imagine all the particles in 
Boutillier's model moving in time in such a way that at each fixed time the picture 
looks like the original Bead kernel model. 

Theorem 1.3. Let a be a real number on the interval (—1,1). In the bulk scal- 
ing limit around a\/2N the Dyson Brownian minor kernel converges to a time 
dependent Bead kernel with parameter a. More precisely, 

(16) irf'="'^((n,x,t),(n',x',t')) = lim e-^(*'-*)(4Ar)^("-"')(2Ar)-5 x 
for 

(17) K^^-'^\{n,x,t),{n\x\t')):= 

- (t)^a"'^{{n, X, t), (n', X', t')) + / ^"'-«e5(*'-*)(«'-2a«)+«(a;-s;) 

where 

(18) 0f^<^((n,x,t),(n',a;',O) = 
'2'^^^-'-') j^H---'{x-y)p^^,_,){y,x' -a{t' -t))dy if {n,t) < {n',t') 
otherwise. 

The limit holds uniformly on compact sets and the contour of integration in { \17\i 
is the straight line between U- and u+ where 



(19) u± = a±iy/l^ 

are two points on the unit circle. 



2 



(The topology used to define compact sets on (N x M^)^ is the product topology 
of discrete topology on N and Euclidean topology on M.) 

Remark that f lTSjl contains the transition density p of Brownian motion, defined 
in f l53|) . rather than p*, defined in fl50|l . This Theorem is proved in Section [HI 
Remark too that, since the K^^^ kernel could only be used along spacelike paths, 
the same is true for our time dependent Bead kernel. Finally notice that the kernel 
^Bead coincides with Boutillier's kernel in |Bou09j with parameter a in the special 
case t = t'. 
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Of course convergence of the kernel in this case does not imply convergence of 
the processes, since the kernel only says something about the behaviour on space- 
like paths. Nor does this Theorem give any hint as to how one might construct 
such a dynamical version of a Bead process. However, specialising the Theorem 
above to t = t' leads to the following Corollary. Though by no means unexpected, 
this result has to our knowledge not previously appeared in the litterature. 

Corollary 1.4. The GUE Minor process, defined in |JN06j . converges in the same 
hulk scaling limit to Boutillier's bead process, defined in |Bou09j . 

Proof. The kernel K^^^ specialised to t = t' is exactly the GUE Minor kernel 
and K^^^'^ with t = t' is Boutillier's kernel. Uniform convergence on compacts for 
the kernels is necessary for process convergence. □ 

The plan of the paper is to first, in Sections [H |3] and SJ outline the necessary 
basic theory about point processes. This is mostly a verbose summary of |BR05t 
IBFPSOTj . Section [6] is devoted to computing certain convolution equalities and 



setting up clever notation so that Theorems II . II and II .21 can be proved at the same 

time, performing the computation only once. In Section [7] the actual computation 
is performed and the article is rounded of by the asymptotic analysis in Section |8l 
The fact that the GUE minor kernel kan be extended to a dynamic version 
in this way begs the question whether something similar can be done with the 
Anti-symmetric GUE minor kernel from |Def08at IDefOSbt IFNOQj . In recent works, 
see lBFP+091 IKSOQj , random walks conditioned to stay in Weil chambers of the 
form 

(20) < Xi < X2 < ■ ■ ■ < Xn 



and their diffusion limits have been analysed. It is reasonable to believe that such 
processes could be realised by Ornstein-Uhlenbeck dynamics on Anti-symmetric 
purely imaginary matrices of odd respectively even size. If so then it appears this 
model could be analysed with the same tools used in this paper and would lead 
to a Theorem similar to 11.11 but with and Anti-symmetric Dyson Brownian minor 
kernel 




and 



(21) 



Xi\ < X2 < ■ ■ ■ < X, 



n 



(22) K- 



■A-VBM 



{{n,x,t), {n,x',t')) 
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That is beyond the scope of this paper, but this kernel is in [BFSOQj recovered as a 
scahng hmit in a certain model related to the totally asymmetric simple exclusion 
process (TASEP). 

While this article was being prepared it came to the attention of the authors that 
Patrik Ferrari and Rene Frings |FF10] were working on related problems. They 
prove an analog of 11.11 for matrices whose elements evolve as Brownian motions 
and also for the Laguerre ensemble. 

2. Point processes 

Let A be a complete, separable metric space with some reference measure A. 
Say M with the Lebesgue measure or Z with counting measure. Let M(A) be the 
set of integer valued and locally finite measures on A. A point process X on A is a 
measure on M(A). It is beyond the scope of this article to give a complete overview 
of the theory of point processes, but some results which we use are detailed here. 

A point process can be represented as 

(23) X = Y,S^. 

where (xj)jg/ are random variables which we shall refer to as the points or the 
particles of X. Think of this as a random configuration of points or particles on 
the space A. In this paper we shall only consider point processes which are simple, 
i.e. all distinct. 

To work with point processes it is convenient to define the so called correlation 
functions. For n = 1,2,..., these are functions p„ : A" — )■ M. When A = Z and A 
is counting measure then 

p„(xi, . . . , Xn) = P[There is a particle at each position Xi, for i = 1, . . . , n]. 

When A = M and A is Lebesgue measure then 

, P[There is a particle in each of \xi, Xi + e), for i = 1, . . . , n] 
Pn{xi, ...,Xn) = lim . 

More generally, see |Joh06] . one can define correlation functions by saying that for 
simple, measurable functions of bounded support, the point process satisfies 

oo „ n 

(24) E[l[{l + <p{x,))] = l + Y,- / l[<p{yj)pn{yi,...,yn)dX''{y). 

i&I n=l "'^ j=l 



A determinantal point process is a point process whose correlation functions 
have the special form 

(25) pn{yi, ...,yn) = det[K{yi, %)]"j=i, 

for some function K : A"^ ^ C This is a very special and simple situation 
since all information about the point process is encoded in this function K of two 
variables which is called the correlation kernel. Nonetheless processes of this kind 
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are commonplace in mathematics today arising from such diverse sources as tihngs 
with rhombuses or dominoes of regions in the plane, random walks, eigenvalues of 
unitarily invariant random matrices and, as shown in the next section, so called 
L-ensembles. Indeed, the main theorems of this paper. Theorems 11.11 and ll.2[ state 
that certain processes are indeed determinantal point processes. Also note that 
with correlations functions of this form, the right hand side of of (l24l) turns out to 
be the definition of the Fredholm determinant of the integral operator with kernel 
K{x,y)(f){y), again see |Joh06] . 

Asymptotic analysis of such point processes can be performed by working with 
the kernels only. 

Proposition 2.1 (Proposition 2.1 in [JN06]). Let X\ X^, , ...be a 

sequence of determinantal point processes, and let X^ have correlation kernel 
satisfying 

(1) —> K , N ^ oo pointwise, for some function K , 

(2) the are uniformly hounded on compact sets in A^ and 

(3) For C compact, there exists some number n = n{C) such that 

det[ii'^(xj,Xj)]i<ij<m = 

if m> n. 

Then there exists some determinantal point process X with correlation kernel K 
such that X^ — )■ X weakly, N ^ oo. 

3. Introduction to L-ensembles 

This section summarises the exposition in [BROS] . The reader who wishes to 
pursue the subject of determinantal point processes will find |Bor09t IMac75] illu- 
minating. 

3.1. Measure theory. Let P, J-") be a discrete probability space. That is, Q 
is a finite set, J-" = 2^ is the cr-algebra of subsets of Q. Furthermore, P : J-" — ^ M 
is a measure satisfying 

(1) F[n] = 1, 

(2) P[0] = 0, 

(3) El, ^2 G J" and n ^2 = implies F[Ei U E2] = P[^i] + P[^2]- 

The last property is called additivity. 

An element G J-" is called an event and F[E] is called the probability of E. To 
sample the distribution P means picking an element u & Q. Since we are working 
with a finite space, ¥[E] can be decomposed as a sum over singleton sets, 

p[E] = 5^p[M]. 
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3.2. Point processes. We now specialise and consider probability spaces of the 
following form. Take a finite set A'. A point process on A" is a probability 
space (f2 = 2'^, P, = 2^ = 2^ ). To sample the point process means to pick 
an element of 2"^, i.e., a subset of X. For compatibility with |BR05j . we will use 



uppercase letters at the end of the alphabet, X, Y, 
elements of Q. 

Given an X G fi, let E;y{X) be the event 



when speaking about 



J^3E;,{X)= U {Y}. 



XCYCX 

The probabilistic interpretation is that Ex{X) is the event that all x ^ X are 
in the chosen set. Note that '■ ^ ^ J-"- 

3.3. L-ensembles. Let us specialise even more. Take a matrix L of size \X\ x 
We shall index the rows and columns of this matrix by X. The L-ensemble on X 
is a point process (fi = 2'^, P, J-" = 2^ ) such that, for X C A", 

det Lx 



(26) 



P[{X}] 



det(l + 

Here Lx means pick out those rows and columns that correspond to X, this giving 
a |X| X |X| matrix. Also, 1 is the identity matrix of appropriate size. For sets in 
J-" which are not singletons, the measure P is defined by the additivity property. 

This is only a well defined probability if the expression f l2B]) is always positive, 
for example if L is positive definite. The fact that the probabilities sum up to one 
is guaranteed by the following well known formula, the Fredholm expansion of a 
determinant. 

Lemma 3.4. Let M be a matrix whose rows and columns are indexed by the finite 
set X . Then 

det(l + M) = ^ detMx, 
xcx 

with the understanding that the determinant of the empty matrix is 1 . 

We do not show this, it can for example be done by induction over the size of the 
matrix. To convince the reader of the validity and triviality of the above lemma, 
let us see what happens with a 2 x 2 matrix. 



det 



1 + a h 
c l + d 



det 



= det 

1 O' 
1 



1 h 
l + d 



det 



1 h 
d 



+ det 
+ det 



h 

l + d 

a 
c 1 



+ det 

det [d] + det [a] + det 
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The first equality comes from tlie fact tliat tlie determinant is linear in tlie first 
and second column. The second equality comes from expanding along rows. 

Theorem 3.5. Let K = L{1 + L)-\ Then for all X G 2^, 

nE;:,{X)] = detKx. 
The matrix K is frequently called the correlation kernel of P. 

Proof. Let be the identity matrix with the ones corresponding to elements in 
X := A" \ X set to zero. 

XCYCX 

= det(l + L)"^ detLy 

XCYCX 

= det{l(^x) + L){1 + L)-' 
= det(l + L-l(x))(l + i:)-i 



det(l-l(x)(l + i.) 



-1^ 



= det[l-(l + L)-i]x 
= det[L(l + L)-i]x 

The third equality is due to a variant of Lemma 13.41 □ 



3.6. Projection on a subspace. Given an L-ensemble on X, let us take an 
arbitrary fixed subset Af C X. Define it's conjugate Af := X \ M . We want to 
study a certain projection P* of P to 2^ . This shall give us a smaller point process 
{n* = 2^,P*, J^* = 22^), specified by 

Again P* is defined by additivity for events that aren't singletons. This should be 
thought of as a conditional probability. Compute for example 

F*[E^{D)] = F*[Udcfcm{F}] 

DCFCAf 
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Theorem 3.7. Let K* = - + L)-'^\^f. Then for all X E 7*, 

Here, means pick out those rows and colums that correspond to M . For proof 
see jBRllH] . 

3.8. Eynard-Mehta Theorem. Again let us speciahse. We are interested in 
studying a point process on the space M = X'^^^ VA X'^^^ VA ■ ■ ■ VA X^^^ which is the 
disjoint union of N finite sets for n = 1, 2, . . . , A^. 
A sample x E M can be written 

™ _ ^(1) 

where a;^"-* G 2'^'"' . Fix an integer p and the following functions: 

0fc M, 



(27) 



for n = 0, . . . , A^ — 1 and = 1, . . . , p. The measure we are interested in has the 
form 



(28) r[{x}] = ^-Met[</>,(a:!°))]f,^,detro(4°\^P0]M=i 



(iV-l) ^(iV)MP 1 ,r ,, / (iV)MP 



■■■X det[iy^_i(xr-^xrOE,.=idet[^,(xrOE,.=i 

if Ix*-"*"^! = = ■ ■ ■ = Ix*-^-*! = p and P*[{x}] = otherwise. Here x^"^^ is the Z:th 
element of the set x^"-* for n = 1, . . . , A^ and / = 1, . . . , p. For that to be well 
defined one needs a total ordering on X^'^\ but which one we use is not important. 
Note that the measure P* : 2^ — )■ M is defined by f l28p for singleton sets and by 
additivity for all other sets in 2^ . Measures of this form turn up everywhere in 
random matrix theory and related combinatorial models. 

The idea now is to extend the space on which the point process lives in such 
a way that the new bigger process admits an L-ensemble representation. Let 
A/" = {1,2, ... ,p} and let X = N' U M . The measure P* can then be expressed as 

^[{xuM}] 



¥*[{x]' 
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where P is the measure defined by fl2^ where 



(29) 



"0 


$ 














-Wo 














-Wi 



































■ 








Here $, Wi, . . . , Wn-i and \l/ are certain blocks and means the zero matrix of 
appropriate dimension. The minus signs are for convenience later. Recall that the 
L-matrix should be of size \X\ x \X\ and that it's rows and columns are indexed 
by elements ofX = AfU X^^^ U X^^^ U ■ ■ ■ U The determinants of the various 

blocks in (l29l) will be exactly the determinants that occur in (l28l) . 

Here, $ and ^ are matrices of dimension p x and lA*^^)! x p respectively 

defined by 

for n e Af and x G 
for X e A'(^) and n eJ\f. 



The matrices Wn for n 
(ITTl)- Let 



1, 



are of size x and defined by 



(30) 



n,m) 



m— 1 ) 



n < m, 
n > m. 



Theorem 3.9 (Eynard-Mehta Theorem). Assume 

(31) M := Wo---W^7V-i^ 

is invertible. Then there exists a correlation kernel K* for the measure P*, that is, 

(32) nEAX)] 
This matrix can be written in block form 



detK*^. 



K* 



K* 

-"■0,0 



K* 



K, 



0,N 



K* 



where 

(33) = Wln,N)'^M-'^Wlo,m) - Wln,m}. 

Note that the block K* „^ is of size l^f^")] x lAf^™)]. 
Proof. We shall use the following matrix identity. 

(34) 



'A 


B 


-1 






C 


D 




D-^CM-^ 


D-^ - D-^CM-^BD-^ 
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Here, A and B must be square blocks and M = BD^^C — A. This is easy to verify 
by explicit computation. 

Now according to Theorem 13.71 we need to invert + L which can be decom- 
posed as the left hand side of ([34D with A = 0, i? = [$, 0, 0, . . . ], 





"1 


-Wo 





. 0" 


-1 p 







1 


-Wi . 


. 




(35) D-^ = 








1 


. 















1 





W[o,i) 
1 




W^[0,2) 
W[l,2) 

1 







W[Q^N) 



'2,N) 



and C appropriately chosen. 
Thus 



BD-^ = [$ <l>l^[o,i) 



[0,JV)J 



D-^C 



and M is as given by (PT]) . Applying (1341) to the formula in Theorem 13.71 gives 

K* = 1-D-^ + D-^CM-^BD-\ 
Inserting the various ingredients above into this formula proves the theorem. □ 

4. Tips and Tricks 

To compute the kernel in the main Theorem, the following additional ideas are 
needed. None of these are new but for the instruction of the reader they are 
summarised here. 

4.1. Continuous state space. In the version of the Eynard-Mehta theorem 
above. Theorem 13. 9[ the state spaces A'("),forn = 1, N, are finite. In the lit- 
erature a version with Af^") = M for all n is more common. We shall now expound 
on the relationship between these two versions of the same useful theorem. 

Again we are faced with analysing a measure on the form fl28|) but now 
P*[{x*}] is the probability density of configuration x*, 
0„, : M — )■ M, for n = 1, . . . , p, 
iy„ : M X M -> M, forra = 1, . . . , AT, 
: M — )■ M, for = 1, . . . , p. 

The state space is now 

(36) MURU---UR = A^xM. 

Pick some discretisation Ai of the real line, i.e. some sequence A^i, Ai2, 
of \M\ real numbers. Restricting the measure in 
fl36|) — to the state space 

(37) MUMU---UM = N X M 



-with state space given by 
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gives a measure on a discrete set of exactly the type to which Theorem 13.91 apphes. 
Now all the blocks in for example fl5S]) are x and thus D is a matrix of 
size (A^IA^I) X (A^IA^I). The idea is of course to take the limit \M\ — > oo. 

The correlation kernel (133|) is computed through suitable matrix multiplications 
and inversions. Consider for example the matrix multiplication of Wi and W2- 

T^iW.W^U. = ^ $^[l^i].,,[W^2],,. ^ / W^ix,y)W2iy,z)dy, \M\ oo 

The constant |A^|^^ can be absorbed into the normalisation constant Z. Thus we 
see that all the matrix multiplications in the expression fl55]) turn into convolutions 
of the corresponding functions in the continuous setting. In later sections we will 
blur the line between discrete and continuous by sometimes using the notation of 
matrix multiplication for convolutions. It is understood that one needs to check 
convergence and integrability when one goes from the discrete to the continuous. 
That offers no problem in our examples so nothing further will be said on that 
score. 

4.2. Unequal number of particles on each level. We shall illustrate this by 
the example of the GUE Minor process. Consider a GUE random matrix M of 
size N X N. Let the eigenvalues of the n x n minor, that is [Mj^]^^-^^, be denoted 

X^^^ > ■■■ > X^n\ Then these vectors A«, A^^) can be seen as random 
variables. 

For this process, the probability measure for all the variables A = (A*^^-*, . . . , A^^-*) 
is absolutely continuous with respect to the Lebesgue measure, thus it has a prob- 
ability density function (p.d.f.). It can be written |BarOH |JN06[ IFN08] as 

(38) p{X) = ^ 1{A« ^ A(2)} ■ ■ ■ l{A(^-i) ^ A(^)} A(A(^)) J] e-^"""'^'- 

n=l 

Here, A^") ^ A("+i) means A^") and A("+i) interlace, i.e. X'^^^^ < A^"^ < A^"+^^ < 
■ ■ ■ ^ '^i'+i^''- We use the increasingly common notation that A denotes the Van- 
dermonde determinant. It turns out that it is practical to introduce fictitious (or 
virtual) variables A^^"* = X^2^ = X^^^ = ■ ■ ■ = X^^ = — oo. Then the interlacing 
condition can be written in terms of determinants |War07] using the Heaviside 
function H{x) = l{x > 0} and the above p.d.f. becomes 

(39) p{X) = ^ det[/7(A?) - Xf%^^,) det[H{Xf^ - Xf%^^,) ■ ■ ■ 

N 

. ■■detmxf^ - Af -))])-., A(A(-)) ^e-(^"'^'^^ 

n=l 

The reader must agree that this vaguely resembles fl28l) except that the dimen- 
sion of the matrices change. The first is a 1 x 1 determinant and the last of 
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size N X N for example. Notice too that the last column of all these matrices is 
identically one because of our choice of fictitious particles above. 

The way to deal with this, first discovered in |BFPS07] , is to form an L- matrix 
similar to the one in fl29l) but which looks like this. 





' 


$ 








... 












-Wo 





... 





L = 


El 








-Wi 


... 







En -I 











... 


-Wn 



Here, in the example of the GUE Minor process, [Wn]x,y = H{y — x) for n = 0, 
1, . . . , — 1, and x, y E M.. For the matrices in f l39|) the last column — which is 
identically one — is moved out to the first column of blocks. Thus a sequence of 
X matrices Eq, . . . , E^-i are produced such that, for m = 1, . . . , A^, 



[En 



n + 1 = m, 
n + 1 ^ m. 



By cranking this machinery it is possible to compute the correlation kernel for the 
GUE Minor process. This is done in |BFPSn7[ IFNOS] . 

4.3. Column operations on the kernel. Recall the expression for the kernel in 
fl33|l . It is sometimes favourable to perform some column operations on the ma- 
trices Win^N)"^ for n = 0, . . . , — 1. Doing column operations means multiplying 
from the right with an upper triangular or lower triangular matrix, say Rm, which 
must be invertible. Typically, but not necessarily, it will have ones on the diagonal 
and a single off-diagonal entry. The kernel ( 133|) then takes the form 

Km = Wln,N)'^RniMRn)-'^Wlo,m) ' W^n,m) 

for 77,, m G {0, . . . , A^ — 1}. 



5. MARKOVNESS ALONG SPACE-LIKE PATHS 

It is a sad fact of life that no transition density for the Warren process is 
known explicitly. We do however know certain marginals. To fix notation let 



X 



[X 



(1) 



,x 



(")) where x'^^^ 



Let I^"-\x) = 1 if the interlacing x*^"^ >-■•■>- 
x*-^-* holds and otherwise. Let Pj:^^ be the transition density for Warren's process 
which we know exists since it is a well defined stochastic process. In this notation, 
|War07| Proposition 6] can be restated as follows. 



Proposition 5.1. For fixed n, x'-"-' and y 

(40) 



in)] 
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where 



e 



(41) pt\x,y) = ^det 

Also, by the characterization after (30) in Warren's paper, it is clear that 

Proposition 5.2. For fixed n, x'^^\ . . . , x'^'^~^\ y^^\ . . . , y^"'~^\ 
(42) 

P(")(x«, . . . , y('\ y(-^) dy^^^ = Pt'^ix^\ x^^'^- . . . , y^^'^). 



We shall compute the eigenvalue measure along a particular space-like path and 
the reader will see how to generalise this. Suppose we want to look at the path 
{n, ti), (n, ^2), — ^2), {n — 2, ^2), {n—2, t^) for some fixed n and < ti < ^2 < ^3- 
The density of the event that the Warren process, started at the origin, takes values 
X, y and z respectively at times ti, ^2 and respectively is by |War07j 

(43) A\xi-^)e'^^M^'^'^^^Pi-^^ 

where t := t2 — ti and t' := t^ — t2- To find the distribution on the afforementioned 
path we need to integrate out x^^\ . . . , y^^\ . . . , y^^~^\ z^"'~^^ and z^'^\ 

We start by integrating out the unwanted x-variables which can be done by 
applying Proposition 15. II which gives 

(44) A^(xW)e-^(^''"')Vl"^(x(«),,W)|j^i^(")(y^ 

Observing that 

l(-\y) = y y y("-^)}/("-2)(y«, . . . ,y("-2)) 



and applying Proposition 15.21 twice we see that 
I{y)Pir\y,z) rfz^")^^^"-^) 

= y y y("-2)}p/r'^(y(^), . . . , z^'\ . . . , z^"-'^) 

We insert that into dH]) integrated and then apply Proposition 15. II to integrate 
out the unwanted ?/- variables to get 

(45) A(x("))e-S(^'"')^||^pS")(xW,yW)l{yW y y y("-2)}x 

A(^ („_2). („_2) j^.2y /--Hz^^\...,z^--^^) 
A(z("-2))^*' ' ' A(;2{"-2)) 

It is well known that 

l{yin) ^ yin-1) ^ det[l{yf) > y'p-'^m,^. 
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if you adopt the convention that ?/n"~^^ = — oo. Then (H5l) can be written as a 
product of determinants. 



(46) A(x("))e-S(^'^"^)^det er'^r' -^^^ 1^ det[l{?/l"^ > yf"^)}]" x 



xdet[l{yr^)>yf"^)}]i:7i,det 



Jm=1 A(z("-2)) 



The same idea can be apphed to any other space-hke path. 

An argument for the corresponding statement for the Dyson Brownian minor 
process is given in |FF10t Section 4]. 

6. Definitions and computations 

Given the theory presented in the last two sections, computing the kernel is 
nothing but a long tedious computation. It was hard to write, hopefully it isn't 
too hard to read. The computation for the Warren process and the Dyson BM 
process can be done at the same time with judicious choice of notation. Let 



(47) i/"(x) : = 



(n - l)!"ix"~il{a; > 0}, n = 1, 2, 
b{x), n = 0. 



be the nth anti-derivative of the Dirac delta function and H := be the Heaviside 
function. 

6.1. Dyson BM. Define the normalised Hermite polynomials, 

(48) hl{x) = [V^n\2'']-'/\-l)''{w*{x))-^D"w*ix) 
which are orthonormal with respect to the weight 

(49) w*{x) = e-''\ 

The transition density of the Ornstein-Uhlenbeck process is the well known 
expression 



(50) pl{x,y) 



7^(1 - qf) 



where 

-t 



While we are at it, define rj" := q* and a*{t) = 1/a/2 for t > 0, and set = (g^ )". 
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6.2. Warren process. For t > 0, let 

(51) h^^\x) = [x/^n!2"r"]-^/2^-l)"(wW(x))-^D"wW(^)- 
These are orthogonal with respect to the weight 

(52) w^'\x) = e-^'/* 

and related to the Hermite polynomials above by h^*\x) — h*{x/\/t). 
The transition density of Brownian motion with variance t/2\s 

(53) p,(a;,^) = ^e-(^--)'/*. 
Define 

(54) = 



s + t 



We are going to set = 1 and a{t) = \/t/2. 

We need the following convolutions in our computations later. 

Lemma 6.3. For n — 1, 2, . . . , and t > 

(55) f h^^\x)w^*\x)p,{x, y) dx = qi'^'''h'^+'\y)w<^'^^\y), 

(56) / h^Ux)w^'\x)H{x -y)dx^ a{t){n + l)-'/'h^\y)w^'\y\ 

JR 

(57) [ H''{x-y)H{y-z)dy^H''+\x-z) 

(58) / pt{x,y)psiy, z) dy = pt+six, z) 
Jr. 

(59) / ptix, y)H^{y - z) dy ^ < / H^{x - y)pt{y, z) dy 
Jr Jr 

(60) / hl^\x)hl^^{x)w^'\x) 
Jr 

(61) pt{x + y,z) = pt{x,z - ny), 

(62) ifA^'i. = itX 

(63) TtTs = n+s 



(64) 



Qt-s n- 



a{s) a{t) 



All of this is also true for the stared functions. The coefficient of x^ in h^n is 



(7"(i)V^!V7r 
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and that is true for h*^ with o replaced by a* . 

These are proved by explicit elementary computation. An important point here 
is that replacing h, q, r, w, p and a with the stared versions these equations still 
hold. By this intelligent choice of notation we can do the computation concerning 
the Ornstein-Uhlenbeck and the Warren process at the same time. Furthermore, 

Lemma 6.4. 

(65) 

r 71—1 7 (s) / w /j-\\k r 

/ H-{x - y)pUy. z)dy = Y: ' L'^'.^ / w^^\y)H-\y - z) dy 

JR ^ rJi^(T(s)v2A;!7r4 Jr 



a(t))" ^ /(/c-n)! x,(t) 



Proof. By orthogonality 

(66) f{x,z) = j H''{x-y)pt-s{y,z)dy 
can be written as 

oo 

(67) f{x,z) = Y,Ck{z)h^,'^ 



fc=0 

for suitable coefficients 



(68) {V2a{s))ck{z) = f{x,z)hi'\x)w^'\x)dx = 

Jr 

(69) = // H''ix-y)pt-siy,z)hi'\x)w^'\x)dxdy. 
Let's start with the case k > n. Apply (159|) . 

(TO) =^ [[ hi'\x)w^'\x)pt.six,y)H^{y - z) dxdy 



Apply ([55]). 

(71) = ^ / h'i\y)w^'\y)H-iy - z) dy 



< t-s 

Apply and (|H7|) n times to get 



(72) = ^(a(«))y^^A«„(^)"-"'(^) 
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Now suppose instead that k < n. Then everything up to flTTl) works the same. 
Apply dsn]) k times. 

(73) = / h^^\y)w^%)H--'{y - z) dy 

rt_s y/kl Jr 

Remember that ho^y) = n^^. 

(74) = / w^Hy)H-Hy - z) dy 

□ 

6.5. Integral representations. With the normalisations above the classical in- 
tegral representations for the Hermite polynomials are 

(75) /i«(a;) = — / v'^e^"-^-^/^' dv 

and 
(76) 



7^4 27ri 



with contours of integration as in Figure [H The starred Hermite polynomials are 
given, for n = 0, 1, . . . , by /i* = hn\ 

As a sort of generalisation of (1751) it can be shown [Nor 09] that for n = 1, 2, 
. . . , 

(77) - 1 = ^ - 
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n = no = ni = S 
^2 = Tt^ = Tt4 = = uq = rij = 2 
"8 = "9 = "10 = "11 = 1 



"12 = 



1* 



9 // 



<3 



Figure 2. Times and levels must follow a space-like path. This 
means that this curve must not take steps upward. 



■ < tN~i and 



7. The Kernel 

The setup now is the following. Pick times and levels (to, '"-o), 
following a space-like path. This means that < to ^ ^ 
no > rii > ■ ■ ■ > njv_i. Without loss of generality we can take njv_i = 1 and 
nm — nm+i G {0, 1} for all m = 0, . . . , N — 2. That is, we end at level 1 and 
only drop one level at a time. For the sake of notation let = and ti\f = t/vr-i. 
Denote by x*-*'-' = {x^i^ > x^2^ > ■ ■ ■ > Xnj!) the eigenvalues at time tk and level 
fik- We shall say that m G Space if the mth step is a down step, i.e. n^, = n^+i + l. 
For fc = 1, . . . , no, let and Sk be the position and time, respectively, of the fcth 
down step. Thus G Space by definition and n,-^ = n + 1 — k. The time of the 
fcth down step is Sk = t^^. = Sno+i-n^^- Also let tq = and sq = to- If the mth 
step is a time step, i.e. Um = Um+i, then we shall write m G Time. That's a lot of 
notation, hopefully Figure [2] should make this clearer. 

Let 



(7J 



(79) WUx,y) 



Ptm+i-t^{x,y) mGTime 
H{x — y) m G Space 



for A; = 1, . . . , no 

for m, = 0, . . . , iV — 1. 



A full configuration of eigenvalues is the A^-tuple 



X 



.x' 



for m = 0, . . . , — 1. We adopt the notation that 



(m+l) 



U 



where x'-™-' G I 

for some large negative real number u if m G Space. In particular = u since, 
by definition iV — 1 G Space. Those are the positions where we step down a level, 
that is, we lose an eigenvalue. One way to think of this is that that particle jumps 
away to some position u which is close to — oo. The weight or probability density 
of configuration x is then, by the Markovness along space-like paths discussed in 
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N-l 



^0) 



MX 



(0) 



(m) ^{'n+l)\-|7i, 



, Xi 



)] 



fc,Z=l' 



m=0 



where we adopt the notation that a^i^+V = u for some large negative real number 



N 



u if Um — Um+i = I, for m = 1, . . . , N. Also xi 
Let the magic begin. Define the block matrix 



u. 



^1) 






$ 














Fo 





-Wo 











Fi 
















Fn-2 











. -Wn-2 





Fn-i 














-Wn-1 


Fn 


















The Fn block in the above block matrix has zero rows but let's keep it notation. 
The {Wn)Z= 

Q are defined above. For / — 1, . 



(82) 
(83) 



= (/)„o+i_;(x) 

Hix — u) 



[Fk]x,i 



0, 



. . . , no, 

k = Tno-l+l, 

otherwise. 



The measure in flHU]) . being similar to that in fl25]) . can then be represented as 
in (l26l) with the above L-matrix. Introduce W[k,i) as in (l30l) . By the general 
theory of these L-ensembles, we need to compute 

(84) K = 1- D-^ + D-^CM-^BD-^ 

for some invertible matrix R. Here, D is as in fl5^ . M = BD^^C, 

(85) 5 = [$ ... 0] , 



(86) 



C 



Fo 
Fi 



The scheme things now is to analyse each of the different components of 
namely BD~^, D~^C, (M)~^ and the upper triangular matrix 1 — D~^. 
First of all, 



= [$ <^Wio,i) $IV[o,2) • • • ^Wio,N)] 
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which is an no x \ matrix. We'll call the kth block of this That is, $o •= ^ 
and 

(87) := ^W[o,k), 

for k = 1, . . . , N, and an explicit expression for it will be given in Lemma 17.51 
Next, let's look at 



(88) D- 





-Fo + iy[o,i)Fi + Wio,2)F2 + Wio,3)Fs 


+ ■■■ + W[o^n)Fn' 








Fl + W^[l,2)i^2 + W^[l,3)i^3 + ■ ■ ■ 


+ W[i^N)Fiy 




^1 




Fn 









where, for /c = 0, . . . , A^, 

N 

(89) ^k = Fk+ ^[k,j)^r 

j=k+i 

As mentioned we will have to do column operations on this, which are represented 
by the tiq x no matrix Rk-, see fl33|) . We will choose Rk in such a way that MRk, 
for A; = 0, . . . , — 1, is asymptotically the identity matrix as u — ?■ — oo. Thus for 
= 0, . . . , A^ let 

(90) := ^kRk. 

These will be explicitly computed in Lemmas 17.71 and 17.81 

Block {k, k') of the kernel in fl84p can with this notation be written, remember- 
ing ([35D, as 

no 

(91) [Kk,k']x,y = -Wlk,k') + J2 \-^k]x,i[M-\A'^k']j,y 

"0 

(92) = -W[k,k') + E i'^kUmRknA^k'hy. 

Note that, as we shall see in Lemma 17. 8[ only the columns 1, . . . , Uk of are 
non-zero. So we need only compute rows 1, . . . , rifc of {MRk)~^. As it happens we 
never need to explicitly write down what Rk is but it can in principle be extracted 
from the proof of Lemma 17. 8[ 

Recall from that M = BD-^C, but what we really need is MRk for k = 0, 
N -1. Multiply dHS]) with ([88]) to compute 

N 

(93) MRk = BD-'CRk = $Foi?fc + ^ <l>Wio,j)F,Rk. 

i=i 

For k = that specialises to MRq = $o^o- For > it will later turn out that 
we only need the first Uk columns of MRk- That allows us to remove those terms 
that only contribute to columns n^ + 1 to rio, which by (jHS]) is those that involve 
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Fq, . . . , -Ffc-1- We shall denote by ~ the operation of removing the unnecessary 
columns. 

N 

(94) MRk - Yl "^^mFjRk 

j=k 

N 

(95) =^Wio,k){Fk+ 

j=k+i 

(96) - 
Another way of saying that is 

(97) [MRk]i,j = [ [^k]i,x[^kUjdx. 
for i, j = 1, Uk. 

7.1. Computing 1 — D^^. From fl35|) it is clear that this is a matrix of size 
A^|A^| X A^|A1|. Furthermore, the {k,l) block of this is identically if / < and 
otherwise —W[k^i). Recall from (!30!) the definition of W[k,i)- For notation introduce 
a function : {0, . . . , — 1} — i- M defined by 

(98) 5(fc):=n^.„o+i--o- 

1=1 

Lemma 7.2. 

(99) [W[k,k')U. = (n,_,J"^'-"*-|^£ff"^-"'='(x-y)p,^,_,,(y,z)rf?/ 

for k < k' and [W[k^k')]x,z = otherwise. 

Recall the definition of H in (1471) . 

Proof. This is shown by induction and k' = k + 1 is the basic case. Let's say 
k G Space. (Recall from the first paragraph of Section [7] what that means.) Then 
Sno+i-rik = tk and rik = Uk' + I so the above expression reduces to [W^[fc,fc+i)]a;,z = 
H{x — z) which is correct. Otherwise k G Time. Then all the S and r factors 
in f l^ equals one. The integral evaluates to [W^[fc,fc+i)]a;,z = Ptk+i-tki^^ ^) which is 
correct. 
Of course 

(100) [W[k,k'+l)]x,z = [W[k,k')]x,y[^k']y,zdy. 

Jr 

If k' G Space then apply first which pops out an rty-tf. = fty-to/''^'tk-to and 
then f l57|) . But k' G Space implies that tk' = Sno+i-n,., = ■5no+i_(„j.,^^+i). On the 
other hand if k' G Time, just a single application of f l58|) completes the induction. 

□ 
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Lemma 7.3. 

(101) 

for k' > k and [W[k,k')]x,z = otherwise. 

Proof. Apply Lemma 16.41 and fl62l) - fl64l) to the formula in Lemma I7.2[ □ 

7.4. Computing BD^^. Recall from fl87|) the definition of which is what we 
shall now compute. Multiply fl82|) with f llOip and use the orthogonality f l60|) to get 



Lemma 7.5. 

(102) = r 



n,,-no 5^ /(nfc^ 01 (to), no-/. (V) f^^^./V) 



for I < Uk' ( equivalently k' > r„+i_/ ) and 

(103) = / w^^-\y)H^-^'iy - .) ci, 

/or I > Uk'. 

7.6. Computing D^^C. We do the same at the other end. Recall the definition of 
in (|89|) . Explicit computations, similar to what we did for $a: in the preceding 
section, leads us to conjecture a general expression, which is proved by induction. 

Lemma 7.7. 

n^.„o+-.-*^ j X j^H-'^^'"\x-y)Ps^^,,.,-tAy.n)dy 

for > I (equivalently k < Tn^^i-i) and [^k]x,i = otherwise. 

Proof. To get the hang of it, consider column one of ^'^ for all A; = 0, 1, .... 
Recall fl83|) . Of course Fn is quite tame, it is identically zero giving us = 0. 
By definition = N — 1, so step number N — 1 E Space and ^n-i = -^Af-i is 
non-zero. Remember that -F/v„i is non-zero only in column one and it is the only 
one of Fq, . . . , -F/v-i that is non-zero in that column. So 

Repeated applications of ( 159|) and fl57|) specialises this to fll04p if you keep track 
of the constants. 
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More generally, for = 0, . . . , A^, recall from ([HSD that column I of is only 
non-zero if k = Tn^+i-i. Remember too from fISI?]) that only occurs in ^q, 
. . . , ^k- That accounts for the fact that [^'fc]^,; = for k > Tn^+i-i (equivalently 
Uk < I). 

The induction on k is now done backwards. Fix / G {l,...,no}. We shall 
show that the formula f ll04p is true for this particular value of /. Specialising to 
k = Tno+i-i, equivalently tj, = Sno+i-i, gives us [^k]x,i = roH^{x-u) = Fj,^ Assume 
now that the theorem gives the correct expression for '^k+i- Since A; 7^ + 1 it is 
clear that Fk+i = in column /. Also, W^fcW^[fc+i,fc = W\^k,k So 

(105) [^k]x,l= [WkUy[^k+l]y,ldy 

Jr 

If k E Time then applying first (IS^ and then flB5]) completes the proof. If A; G Space 
then flFr|) does the same. The calculation is in reality only a question of keeping 
track of the coefficients. □ 

We now exercise the right to do column operations. Recall from fl90p the defi- 
nition of "^k- Rk is chosen judiciously to give a very simple expression for MRk- 

Lemma 7.8. There exists an Rk such that for x E Ai and u a large negative real 
number, 



(lObj [^k\.,i-K,.Ax)r,^^,^^q,^,,^ Y(n.-/)!S(0) V2a{tk) ^ ' 

for 1 = 1, Hk, and ['^k]x,i = fori = nk + l, . . . , Uq. 

Proof. This bit is rather hairy actually. Fix k and let pi = Sno+i-n^+i ~ tk- Let 

flix,u) = C{l)[^k\x,n,-l 

for / G {0, . . . , rifc — 1} where C{1) is the constant that gives 
fi{x,u)= [ H^^\x -y)ppXy,u)dy 

H^+^{x - y)pp^ iy - r-p^u, 0) dy 



/ H^+^{x~y-rp^u)ppXy,Q)dy. 



The second equality is due to fl6Tl) and the third one is a change of variables 
y ^ y + r~^u. Now what we need to do is to produce a sequence of functions go, 
. . . , grik-i such that for each /, 

(1) gi is a linear combination of /o, . . . , fi, 
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(2) gi{x, u) differs from polynomial in x of order / by at most a term exponen- 
tially decreasing in u. That is, there is a polynomial hi {x) of degree I such 
that 

gi{x,u)^h{x)+0{e-''"'''). 
Just to see how it works consider 

(107) /o(a:, u) = I H{x -y- r^^u)pp^{y, 0) dy 

(108) = J pp,{y,0)dy + 0{e-^'/') 

as w — >■ — oo. So let giQ = /o implying that ho is a constant. 
Now for some constant 

" (/ ^p^^^'^^^y) (/ ppo^y^^^^y^ 

look at 

fi{x, u) - Cr~^^ufo{x, u)^ j {x-y- r~^u)H{x -y- r~^u)pp^ {y, 0) 

-Cr-^K J H{x-y -r-^\)pp^{y,0) 

(where C is chosen so that, in the limit of large negative u, the integrands partially 
cancel out.) 

= j{x-y)pp,{y,0)dy + O{e--'/') 

J Pp,{y,0)dy- J ypp,{y,0)dy + O{e-'''/^). 

Thus set gi {x, u) — fi — Cr^^ufo and hi accordingly. 

Now that we see it works for / — and 1, let's do the induction. Say the 
statement is true for /, then 

i=o ^ ^ 

The first term is already a polynomial in x of degree / + 1 modulo an 0-term. The 
other terms are polynomials in x of lower degree which by the induction hypothesis 
can be expressed in gg, . . . , gi. 

What we have effectively done is that we have taken various linear combinations, 
encoded in R, of columns of producing 

['^kR]x,i — polynomial of degree I + perturbation exponentially small in u. 
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Now we will go ahead and do more column operations and multiply the columns 
with suitable constants. Let B! be the matrix that encodes the column operations 
and multiplications that turn the columns of into Hermite polynomials with 
factors as in the statement of the Lemma. Set Rk = RR' and \E'fc = ^kRk and we 
are done. □ 

7.9. Computing (MRk)^^. First observe that M is essentially lower triangular. 
To see this, for I < i < j < no consider 



(109) 



0(e-"'/2) 



This is a computation. You insert the expressions for $o and "^q, apply first fl57p . 
then then fl5^ and lastly perform riQ—j applications of fISBl) . By (I^Ul) the conclusion 
in (UnnD holds. 

Let /„ be the identity matrix of size n x n. In block form M and Rk are: 



(110) 



MRk 



* 0(e-"'/2) 





'* 






Ino-nk_ 





0(e-"'/2) 



The top- right block of MRk is essentially by construction. To see that, re- 
call fl97j) , insert Lemmas 17.51 and 17.81 and apply fl60l) . To invert this matrix use flMl) 
with M = —Ifik + 0(e~"^/^) and B = 0(e~"^/^) and we are only interested in the 
top two blocks. 



Proposition 7.10. 

(Ill) [{MRk)-%, 
fori = l, . .. , Uk and j = 1, . . . , Uq. 



Essentially, it is diagonal (Hurrah!!) due to our judicious choice of polynomials 
and furthermore the identity matrix due to our choice of factors in Lemma 17.81 



7.11. Putting it all together. Recall the expression for the kernel in For 
the case k > k' , inserting the results of Lemmas 17.51 and I7.8[ (I62|) - (l64|) . taking 
account of the cancellations and letting m — )■ — oo gives 



(112) K{{nk, X, tk), {nk',x', tk>)) = r^^.^^ 



Sjk) iaitk>)r---^' 
S{k') V2a{tk) 



1=1 



{uk-iT'^' 



In this case you don't need to use f llOSp . 

In the case k < k' we also need Lemma [7.31 The first Uk terms in the summand 
in ( I92p cancel the corresponding terms in f llOip and we are left with the following 
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infinite sum. 

(113) ir((nfc,x,tfc),(nfc.,x,tfcO) = -r^^-to ^^y-^ Wafa) ^ 



E\/{^<:'^'"r'e,(-)C;^<(-')"-'''''(:^'). 



l=—oo 

We shall now specialise to the case of Dyson BM and Warren respectively. 

Proof of Theorem Insert the various expressions in section [6?T] into f lll2p and (I113p 
and multiply with the conjugating factors 

^ ' 5(A;)e"(*-*") 

which always cancel out when you take determinants. 

To show that the integral representation equals this shifted sum of Hermite 
polynomials, observe that 

1 V V \ / U 



E 



1=0 

and plug in the integral representations in Section 16.51 The function should be 
expanded according to Lemma 16. 4[ □ 

Proof of Theorem Insert the various expressions in section I^T^ into flll2p and flll3p 
and multiply with 

'^(feQ g(n'-«)/2 
^^^^^ S{k)^ 

which always cancels out when you take determinants. The integral representation 
is done exactly as in the last proof. □ 

8. ASYMPTOTICS 

We shall now study the scaling limit of the Dyson Brownian Minor kernel 
from fllOp in the bulk. 

Proof of Theorem \l.'J\ . This is done by saddle point analysis. Let r^r = e^''*'^*^/^^. 
We start out with the case (n, t) > {n',t'). Use the representation ffTOjl and sub- 
stitute u ^ u and V I— 7- TJ^^v^y N/2 in the integral. We must now compute 

1 r r -r'N-n'-l n' 

(116) lim du dv^ x 

Af^oo 27ri Jy V - u u"- 

X exp [{N/2){t^'^v'^ -u^ + 4a(M - r^^v) + 2lnv - 2lnu) + ux - r^^vy] . 
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The part of the exponent that is multiphed by {N/2) is f{u) — f{v) for 
(117) f{u) = -u"^ + 4:au-2\nu. 

The saddle points, satisfying the equation f'{u) = are given by (fT^ and are both 
on the unit circle. We let 6 be the angle specifying those points, so that u± = e^*^. 
The Taylor expansions around these critical points are 



:ii8) 



f{u± + h) = f{u±) + {T^aVl - a? + 4(1 - a^)i)h^ + o{h^). 



The only thing we really need is that the coefficient of h? is non-zero. 

Now deform the contour F to which is a straight line from a — ioo to a + ioo. 
See Figure [TJ When deforming the f -contour through the w-contour out pops the 
residue at u = v which is exactly 



(119) 



— ^n'-n^\{t'-t)(u^-2au)+u{x-y) 



It remains to show that all the other parts of the contours evaluate to zero. 

To see what happens around the saddle point m+ perform the change of variables 
M (— 7- + u/ \fN and v ^ + v / \/N. 



du / dv 



T 



-N-n'-l 



N 



1 

2ni J J V — u 
"1 
_4 

which, as — 7- oo, tends to 



X 



X exp 



(m+ + u/VNY 
f"{u+){u' - v^) + u+{x - T-'y) + 0(1 - rjv) 



0{N- 



1 

27n 



du / dv 



V — u 



n+)"-"exp 



v^) +u+{x - y) 



Switching u i— )■ —u and v i— )■ —v makes the integrand change sign. Thus the integral 
must be zero. The saddle point at m_ contributes zero by the same argument. 

For the remaining contours we need to show that g{u, v) := Re f{u) — f{y) < 1. 
Well, if g{u,v) > 1 somewhere then the integral would not be convergent. But 
the whole thing is a probability density and must therefore be finite, so this is a 
contradiction. If at some point (m, v) it happens that g{u, v) = 1 then f'{u) ^ for 
we have already accounted for all saddle points. Thus by the definition of derivative 
there must be a point nearby where g{u,v) > 1 which is a contradiction. 

Next we need to compute the same scaling limit of the function from (11 II) 
which is a very straight forward computation, yielding f[T5]) . □ 
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